Scalar dark energy models mimicking ACDM with arbitrary future evolution 
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Dark energy models with various scenarios of evolution are considered from the viewpoint of the 
formalism for the equation of state. It is shown that these models are compatible with current 
astronomical data. Some of the models presented here evolve arbitrarily close to ACDM up to the 
present, but diverge in the future into a number of different possible asymptotic states, including 
asymptotic de-Sitter (pseudo-rip) evolution, little rips with disintegration of bound structures, and 
various forms of finite-time future singularities. Therefore it is impossible from observational data 
to determine whether the universe will end in a future singularity or not. We demonstrate that the 
models under consideration are stable for a long period of time (billions of years) before entering a 
Little Rip/Pseudo-Rip induced dissolution of bound structures or before entering a soft finite-time 
^"t ■ future singularity. Finally, the physical consequences of Little Rip, Type II, III and Big Crush 

singularities are briefly compared. 
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I. INTRODUCTION 



The discovery of the accelerating expansion of the universe [l|, l2| has raised a number of difficult problems in 

55)' cosmology. The cosmic acceleration can be explained via the introduction of so-called dark energy (for recent review, 

see |3|,y|) with quite strange properties like negative pressure and/or negative entropy, invisibility in the early universe, 

etc. According to the latest supernovae observations, the dark energy currently accounts for 73% of the total mass 

energy of the universe (see, for instance, Ref. j5(). 

iq , The equation of state (EoS) parameter wu for dark energy is negative: 

Q\ ; wr> =Pd/pd < 0, (1) 



where p^> is the dark energy density and pu is the pressure. Although astrophysical observations favor the standard 
ACDM cosmology, the uncertainties in the determination of the EoS dark energy parameter w are still too large to 
p^i define which of the three cases w < -Im = —1, and w > — 1 is realized in our universe: w = — 1.04+q'5q [flQ- 

y—( ' If w < — 1 (phantom dark energy) [8J], we are dealing with the most interesting and least understood theoretical 

case. For a phantom field, the violation of all of the four energy conditions occurs. Although this field is unstable 
from the quantum field theory viewpoint [9|], it could be stable in classical cosmology. Some observations [10[ may 
be understood as the indication of the crossing of the phantom divide in the near past or in the near future. A 
very unpleasant property of phantom dark energy is the Big Rip future singularity [8|, Illl - ll4| , where the scale factor 
becomes infinite at a finite time in the future. A less dangerous future singularity caused by phantom or quintessence 
dark energy is the sudden (Type II) singularity [15[ where the scale factor is finite at Rip time. Nevertheless, the 
condition w < —1 is not sufficient for a singularity to occur. Mild phantom models where w asymptotically tends to 
— 1 and the energy density increases with time or remains constant but there is no finite-time future singularity are 
discussed in recent works fla - [l8j . The key point is that if w approaches —1 sufficiently rapidly, then it is possible to 
have a model in which the time required for the occurrence of the singularity is infinite, i.e., the singularity effectively 
does not occur. However, if the energy density grows, the disintegration of bound structures eventually occurs in a 
way similar to the case of the Big Rip singularity. Such phenomena, called the Little Rip and the Pseudo-Rip, destroy 
the bound structures in the universe at a finite time! 

The most convenient formalism to construct the (scalar or fluid) dark energy is to use the EoS: 

P = 9(P) , (2) 

where g is a function of the energy density. The evolution of the universe then depends on the choice of the EoS. 
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The aim of this article is to develop a general approach for the construction of dark energy models which are 
compatible with observational data and which provide various scenarios of evolution. Moreover, we arc most interested 
in models which are stable for a long period of time before entering the Little rip/Pseudo-Rip induced dissolution 
of bound structures or before entering the mild finite-time future singularity (with a finite scale factor at the Rip). 
This question is analyzed from the viewpoint of the dark energy EoS and corresponding description in terms of 
scalar field theory. In Sec. lU the general approach to this problem is developed. The classification of cosmological 
models with different future evolutions is presented. The confrontation with observational data is also given. The 
non-singular Little Rip cosmology based on a massive scalar potential is considered in Sec. IIIII It is demonstrated 
that it may be indistinguishable from the ACDM cosmology, being stable for a long time before the disintegration 
of bound structures. Although the similar models are considered in [Tg|, it is interesting to note that the simplest 
model with Klein-Gordon potential is described apparently for the first time. In Sec. lIVI the asymptotically de Sitter 
quintessence or phantom cosmology is proposed. The properties of such realistic cosmologies are similar to those of 
the previous section. In comparison, for example, with [17J, our approach is based on equation-of-state formalism 
instead of setting of the Hubble parameter as function of time. The next two sections are devoted to the construction 
of quintessence dark energy with a Type III future singularity and phantom/quintessence dark energy with a Type II 
future singularity. These results based on EoS formalism are novel in cosmology. The estimation of time remaining 
before Type II or Type III singularity is made for first time. The comparison of the predictions of such models with 
observational data demonstrates that they may be indistinguishable the ACDM model up to the present. Moreover, 
they may be stable for billions of years before reaching a future singularity. In this sense, such models represent a 
quite viable alternative to ACDM. Sec. lVIIl is devoted to the construction of Big Crush scalar quintessence models and 
their comparison with the Little Rip or Type II future singularity cosmology. Specifically, some physical properties of 
the Big Crush versus the Rip are discussed. Phantom models are briefly reviewed in Sec. I Villi The main qualitative 
result of our study is that the current observations make it essentially impossible to determine whether or not the 
universe will end in a future singularity. We should note once more that the consistent application of equation-of-state 
formalism is the main feature of our work. Some summary and outlook are given in the Discussion section. 

II. SCALAR DARK ENERGY MODELS 

In this section, we consider the construction of one-scalar dark energy models and compare the models with 
observational data. We mainly concentrate on different phantom theories whose classical behavior is very similar to 
that of the ACDM model. 

For the spatially-flat FRW universe with metric 

ds 2 = dt 2 -a 2 (t)(dx 2 + dy 2 + dz 2 ), (3) 

the cosmological equations, that is, the FRW equations are given by 

° V = f- >--»(!) o + »- <«> 

where p and p are the total energy density and pressure, a is the scale factor, ' = d/dt, and we use the natural system 
of units in which 8ttG = c = 1. 

It is convenient to write the dark energy equation of state (EoS) in the following form: 

Pd = -pn - /(pd) , (5) 

where /(pd) is some function. The case /(pd) > corresponds to the EoS parameter w < — 1 (phantom) while the 
case f(p) < corresponds to the EoS parameter w > — 1. If dark energy dominates, one can neglect the contribution 
of other components (matter, dark matter). Then from Eq. ([4]), one can get the following expression for time variable: 

2 f x dx 

t = T*Lwr x = ^ p - (6) 

Hereafter, it is convenient to omit the subscript D. The quintessence energy density decreases with time (x < Xq), 
while the phantom energy density increases [x > xq). Following Ref. |l4[, one can find the following behavior for the 
expression ((6]): 



1. The integral ([6]) converges at p — >• oo. Therefore a finite-time singularity occurs: the energy density becomes 
infinite at a finite time if. The expression for the scale factor 

/ 2 f x xdx \ 

a = aoCxp Ui7Rj' (7) 

indicates that there are two possibilities: 

(a) The scale factor diverges at a finite time (Big Rip singularity [a. UlL Il3| ) . 

(b) The scale factor remains finite; however, a singularity (p — > oo) occurs. This is a Type III singularity J19J . 

The key difference between ([Taj) and (|lb[) is that for (|lb[) the energy density grows so rapidly with time that the 
scale factor does not reach an infinite value. 

2. The integral ((6]) diverges at p -> oo. Such models are described in [16[ (see also Refs. [18|, |2(|). The energy 
density grows with time but not rapidly enough for the emergence of the Big Rip singularity. According to 
Ref. [16| , we have a so-called "Little Rip" : eventually a dissolution of bound structures occurs at a finite future 
time. Nevertheless, formally the future singularity does not occur (or, rather, it is shifted to the infinite future). 
Such scenarios are possible only in the case of phantom dark energy. The next two scenarios are possible for 
both phantom and quintessence dark energy. 

3. The integral (|6|) diverges at p — > pf. The dark energy density asymptotically tends to a constant value ("effective 
cosmological constant" ) . Such asymptotically de Sitter theories represent the natural alternative to the ACDM 
model, which also leads to non-singular cosmology. Nevertheless, even for a non-singular asymptotically de 
Sitter universe, the possibility of a dramatic rip which may lead to disappearance of bound structures in the 
universe remains possible [ly, [l8[ . 

4. Another interesting case corresponds to f(x) — > ±oo at x = xt , i.e., the pressure of the dark energy becomes 
infinite at a finite energy density. The second derivative of the scale factor diverges while the first derivative 
remains finite. It is interesting to investigate the properties of dark energy with such a (sudden or Type II) 
finite-time future singularity [15l . [l9( . 

What does this last type of singularity mean from the physical viewpoint? As the universe expands, the relative 
acceleration between two points separated by a distance I is given by la/a. If there is a particle with mass m at each 
of the points, an observer at one of the masses will measure an inertial force on the other mass of 

F inor = mla/a = ml (h + hA . (8) 

Let us assume that two particles are bounded by a constant force Fq. If F[ ner is positive and greater than Fq, the two 
particles become unbounded. This is the rip produced by the accelerating expansion. Note that Eq. (|5|) shows that 
the rip always occurs when cither H diverges or H diverges (assuming H > 0). The first case corresponds to the Big 
Rip singularity, while if H is finite, but H diverges with H > 0, we have a Type II or sudden future singularity, which 
also leads to a rip. Even if H or H goes to infinity at the infinite future, the inertial force becomes larger and larger, 
and any bound object is ripped, i.e., the Little Rip cosmology emerges. If H is finite and H is negative and diverges, 
then all the structures are crushed rather than ripped. 
Eq. ([8]) can be rewritten as 

F ineI = ml(^- + ttp) . (9) 

The crush corresponds to f(x) — > — oo at x — > Xf < oo while the case f(x) — > oo at x — > X{ describes the sudden 
future singularity. The equivalent description in terms of scalar theory can be derived using the equations: 

p = ±tf/2 + V(cp), p = ±4> 2 /2-V(0), (10) 

where <ft is the scalar field with potential V(4>). The sign "— " before kinetic term corresponds to the phantom energy. 
For the scalar field and its potential, one can derive the following expressions: 

^'^tJIJwv (11) 

V(x)=x 2 + f(x)/2. (12) 



Combining Eqs. (fTTj) and (fT2|) gives the potential as function of the scalar field. For simplicity, we choose the sign "+" 
in Eq. ([TTj) for phantom energy and "— " for quintessence. For the crush and sudden future singularity the potential 
of the scalar field has a pole, i.e. from the mathematical viewpoint, these singularities are equivalent. Note that 
singularities often correspond to the infinite value of the scalar field 4> — > ±00. For the sudden future singularity 
potential, we find V(4>) — > +00, and for the big crush, V(</>) — > —00. 

Confrontation of the theoretical models with observational data consists mainly of comparison with the distance 
modulus as a function of redshift from the Supernova Cosmology Project [28j . The distance modulus for a supernova 
with redshift z — ao/a — 1 is 



fi(z) = const + 5 log D(z) . 



(13) 



where D{z) is the luminosity distance. As is well-known, the SNe data are well-fit by the ACDM cosmology. For such 
a model (which we call the "standard cosmology" (SC)), we obtain 



£>£ C - 7r(l + z) / (On(l + zf + «a) 
-"0 Jo 



-1/2 



dz . 



(14) 



Here, £l m is the fraction of the total density contributed by matter, and fl\ is the fraction contributed by the vacuum 
energy density. 

One can get also the deceleration parameter go and jerk parameter jo |2l| : 



9o 



Jo = 



1 


d 2 a 


aH 2 dt 2 


1 


d 3 a 



1H 3 dt 3 



t=t 



t=o 




N=0 



1 d 2 (H 2 
2H 2 dN 2 



3 d(H 2 
2H 2 dN 



(15) 
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Here N is defined by 



N = - In (1 + z) 



(16) 



For the current time t = to, we have N 



d x d 

dN = - {l + z) dz 



0. It is useful to note that 

d 2 



dN 2 



, d , ,9 d 2 



(17) 



Measuring the deceleration parameter and especially the jerk parameter is a much more difficult task than measuring 
Hq. In order to measure the Hubble constant, one needs to derive the distances to objects at ~ 100 Mpc; this 
corresponds to a redshift of z > 0.02. To obtain go, one needs to observe objects to redshift z > 1. Therefore, current 
observational results for the deceleration and jerk parameters are not totally reliable. For example, Ref. [22| gives for 



a flat model tight constraints on go = —0.81 ± 0.14 and jo = 2.16_ 'y 6 from type la supernovae and X-ray cluster gas 
mass fraction measurements. These results are consistent with ACDM at about the la confidence level. 
For the ACDM model, one gets 



9o — :^ m — 1 , 



Jo 



1. 



(18) 



Therefore, measuring Jo is important in the search for deviations from ACDM, since all ACDM models, regardless of 
the matter and cosmological constant energy densities, are characterized by j = 1. 

This concludes our general discussion of scalar dark energy models. Below we present several examples of such 
models. 

III. SCALAR LITTLE RIP COSMOLOGY 

It is known that the Little Rip cosmology can be realized in the class of exponential or power-law scalar potentials 
[I3, LL8(- ^ is interesting that even for the case of the simplest Klein-Gordon potential, the Little Rip can occur. 
Let us start from 



f{p) = a 2 = const . 



(19) 



From the conservation law for the dark energy fluid ((4j) , it is easy to obtain 

p = p -3a 2 \n(l + z). (20) 

The luminosity distance is given as 

D L = —(l + z) (n m {l + zf + (l + 3(w + l))n D y 1/2 dz. (21) 

-"o Jo 

where the current EoS parameter wq is 

w = -1 - — , (22) 

Po 

and Q m and Qb are the matter and dark energy contributions to the total energy budget. It is clear that if a 2 <C po 
(and therefore Wv sa —1) there is good agreement with the observational data: such a cosmology is indistinguishable 
from the ACDM model. The deceleration and jerk parameters are given by 

3 3 9 

q = -a m -l+ -(wo + l)n Dl jo = 1 - 2^0 + 1)«d • (23) 

For —1.15 < wo < — 1 and ^d = 0.72, these parameters are in good agreement with observational data. The scalar 
field grows with time as 

= 0o + at , (24) 

and the potential describes the massive scalar field: 

^) = ^(0-0*) 2 + Y, m 2 EE3« 2 /2, 0* = ^_|^. (25) 

This type of potential is very popular in particle physics. For example, the light scalar particles like dilaton and 
moduli appear in superstring theories. The dark energy density increases with time. Hence, the universe accelerates. 
One can estimate the time required for disintegration of the Sun-Earth system, as an example. The dimcnsionlcss 
inertial force 



(26) 
can be expressed for t >• t as follows 



aH$ 



Finer ~ ^D ■ (27) 

Po 



Taking into account that the phantom energy density increases with time as 

Po /2 -^-K + lw) , (28) 

and that the Sun-Earth system disintegrates when F incr ~ 10 23 (see [3), one can find that the disintegration time is 

10 12 10 13 

\wo + l\H |t«o + l| 

Hence, we have presented a realistic Little Rip cosmology caused by scalar dark energy with a standard particle 
physics massive potential. Note that if such cosmology occurs, one can speculate on visible reduction of galaxy 
clusters number in future cosmological surveys. 



IV. ASYMPTOTICALLY DE SITTER EVOLUTION: PSEUDO-RIP 

It is evident that some phantom/quintessence models may describe asymptotically de Sitter evolution. However, 
even in this case the disintegration of bound structures may take place. Such a scenario was dubbed the Pseudo-Rip 
cosmology [ly, [l8|. Of course, not all asymptotically de Sitter phantom scenarios lead to disintegration of bound 
structures. 

Let us consider the example of a fluid which describes the phantom and/or quintessence held asymptotically 
approaching the de Sitter regime. Let 

/(p)-« 2 sin^. (29) 

If 2kpf < p < (2k + l)pi, k = 0, 1, ..., we have a phantom case while for (2k + l)pi < p < (2k + 2)pi Eq. (129)) describes 
quintessence. For p <C pt , the parameter w is nearly constant: 

2 

a it 

W BS — 1 , 

Pi 

P = Pi [ ±- arctan f tan f J J + 2fc J , 6 = ia 1 j-Kpi , A = p /pf , (30) 

where the "±" corresponds to phantom and quintessence theories, respectively. Depending on the parameter A, the 
dark energy density tends to a single value from the set of "effective cosmological constants" 

A cff = (2k + l)p { . (31) 

If the effective cosmological constant is sufficiently large (a /a 3> 0), then disintegration of bound structures can occur. 
It is obvious that this scenario occurs only for the phantom case (for quintessence with asymptotic de Sitter evolution, 
the acceleration of the universe can only decrease) . The dimensionless inertial force 

a 

F iner = <3 ■ V^ / 

p a 



tends to pf/po = A -1 (for k = 0) because a ~ exp(-y/pf/3t) at t — > oo. Therefore if A < 10 -23 then F incr > 10 23 at 
t — > oo and the disintegration of the Sun-Earth system eventually happens. 

Such a scenario is compatible with observational data. For small A, one can write with good accuracy the past 
dark energy density 

p^p (l + z)- s . (33) 

For small 5 (w — > —1), the dark energy density is nearly constant in the observable range < z < 1.5. Therefore, the 
observable relation between modulus and redshift can be fulfilled in the model given by Eq. (|29p . 
The deceleration and jerk parameters are within the range of the standard cosmology (as S — > 0): 

go = -O m --Q D -l, (34) 

io = l + ^^n D . (35) 

Thus, a realistic mild phantom scenario (Pseudo-Rip) may be easily realized. 

V. DARK ENERGY MODELS WITH A TYPE III FUTURE SINGULARITY 

We shall consider in this section a flat universe which ends up in a Type III future singularity [19(. Let us start 
from 

g(p) = -{3 2 a\p l +^ , (36) 

so that 

f(p) = p(-l + p 2 a\p^), 



where (3, at, and e are positive constants. One can find the dependence of the dark energy density on the scale factor 

p = /T 6/e (of - a e )~ 3/e . (37) 

Putting, for example, e = 1, we have 

p=-f3 2 a iP ^. 
In this case, one can find the scale factor in parametric form 

a — <2f sin 77 , 
1 



t = if + - ( In 

K 



tan 



'/ 



1 



where 



1 



2sqrt3f3 



-3/2 
3 a f ' 



cos ?y 



dt 



■ cos 77 , 



(38) 



cos 4 77 
k sin 77 



di] . 



We now set if = 0. Therefore 77 = corresponds to t = —00, 77 = tt/2 to i = (future singularity) and 77 = 7r to 
i = +00. Hence, this solution describes two universes: the first one begins at i = —00 (Big Bang) and then expands 
to a future singularity which takes place at i = 0. The second solution begins at t = (at a singularity) and then 
progressively contracts until a big crunch singularity at t = 00. The asymptotic behavior of the scale factor near the 
future singularity in both cases is 



a = ai (l-(5 K ) 2/5 |i| 2/5 ), t~0. 



The addition of dark matter allows us to construct cosmological models in which the age of the universe is close to 
the conventional value of 10-20 Gyr. The dependence of the dark energy density on redshift is given by 



p = P o{l + zf 



No-1 



:!/«■ 



K N (i + zy-i. 

where jVq = (<2f/ao) e . For the current value of the EoS parameter wq we have 



w 



-(3 2 N a alpf = - 



N a 



No -I 



(39) 



(40) 



One can use the standard relation between redshift and time 

dz 



Ho~ l 



= - dt 



(l + z)^h(z) 

h{z) = n m (i + zf + n D (i + zf 



No -1 



No(i + z y-i 



3/« 



(41) 



for the calculation of the age of the universe and the estimation of the time of the future singularity. Integrating 
Eq. (|4"Tj) from z = 00 (i = 0, Big Bang) to z = (i = i ) gives the age of the universe: 



tfo" 1 



dz 



/O (l + z)y/h(zj 

For No ^$> 1 (i.e., for wo ~ —1) the function h(z) can be approximated by 

h( z ) w n m (i + z) 3 + n D . 



(42) 



(43) 



Therefore, the age of the universe is eventually independent of e. This parameter, however, may change the remaining 
time before the future singularity if — to- Note that for the calculation of this time we can use Eq. (14"T1) . simply 



Table I: A numerical calculation of the age of the universe, to, and the difference between the future singularity time, if , 
to for various values of Wo and e. The time unit is 10 9 years (Gyr) and we choose Hg 1 = 13.6 Gyr. 



and 





Wo = -1.01 


TOO = -1.05 


wo = —1.1 


e 


if — to 


to 


if — to 


to 


tf — t(j 


to 


1 


52.95 


13.66 


30.48 


13.73 


22.17 


13.81 


2 


29.45 


13.66 


17.79 


13.71 


13.31 


13.78 


5 


12.64 


13.65 


7.93 


13.69 


6.08 


13.73 


10 


6.40 


13.65 


4.09 


13.67 


3.17 


13.69 


50 


1.27 


- 


0.83 


- 


0.65 


- 


100 


0.63 


- 


0.41 


- 


0.33 


- 


1000 


0.06 


- 


0.04 


- 


0.03 


- 



assuming that the variable z can take negative values. The lower limit of integration corresponds to the scale factor 



Of, i.e., Z[ = N, 



-l/e 



1. Therefore for if 



to one gets the following relation 
dz 



if — to — H Q 



Zi {l + z)s/h(z) 



(44) 



It is obvious that our model can fit the Supernova Cosmological Project data. For N 3> 1 the dark energy density is 
nearly constant in the interval < i < to, i.e. the model (|36[) mimics a cosmological constant in the past but leads 
to a finite-time future singularity. Moreover, the observational data do not impose any significant restrictions on the 
lifetime of the universe. 

The numerical calculation of the age of the universe, to, and the difference between the future singularity time, if, 
and to for various values of u>o and e are presented in Tabic HI The value of the Hubble parameter is chosen to be 
13.6 Gyr for this calculation. 

The deceleration parameter is found to be 



tfo" 1 



3 3 

qa = 2^™ - 1 + 2^d(1 + wo) 



and for wq ~ — 1 qo ~ 9o • For the J er ^ parameter we have 



jo = r(3 + e )( w o + u>o)^D + l 



(45) 



(46) 



For sufficiently small e and w ~ — 1, the jerk parameter is nearly equal to 1. Hence, a viable quintessence dark energy 
model which is compatible with observational data and leads to a Type III future singularity is constructed. Of 
course, it may be presented in terms of a scalar field with the field value <f> and potential V(<fi) written as functions of 
r\. Furthermore, the dependence of V(<j>) on <j> can be presented explicitly. We have 



V_ 
no 



+ 00) 



V{4>) = ai(cos(aoW> - </>o)))" 6A (l + co S - 2 (a o (0 - <fo))) 



a 



e , 



OL\ 



\P 



-6/e, 



(47) 



Other models with similar properties can be constructed. 



VI. TYPE II FUTURE SINGULARITY DARK ENERGY 

In this section we discuss realistic models of dark energy which contain a Type II future singularity 11911 Tor sudden 
future singularity [15|). It is known that such evolution may be also realized in f(R) modified gravity [23| . 
The simplest choice for an EoS producing a Type II future singularity is 



f(p) 



1 - P/Pi 



(48) 




Figure 1: The scalar potential for the model (|48|l. For a crush the scalar field rolls down from — oo to and V((j>) —¥ oo; for a 
phantom sudden future singularity the scalar field rolls up and V(0) — » oo at some (f> — cj> s . 



where a and p{ are positive constants. For po < pi , such a model describes phantom energy. Its energy density grows 
with time until the pressure tends to infinity and a phantom sudden future singularity occurs. If po > Pi-, the energy 
density decreases and a big crush occurs when p is equal to pf. 
The time remaining before the future singularity is 



if — t 



V3 Jxa 



d.v 



(49) 



The corresponding description in terms of scalar field theory can be derived through Eqs. (fT2|) and (fTT|) . The potential 
of the scalar field in parametric form is: (i) for phantom energy 



0(2/) 

V(y) 



== (arcsiny + y{\ - y 2 ) 1/2 ) 



2y 2 



l-y* 



, < y < 1 , 



(ii) for quintessence 



4>{y) 

V(y) 



v^? 



3n 



2y + V-i 



ln(l - v/l -y 2 ) +lnx/J , 

0<y<l, 



(50) 

(51) 

(52) 
(53) 



and 7 = a 2 j pf. This potential is depicted in Figure [TJ 

Such a model in principle can fit the latest supernova data from the Supernova Cosmology Project. The dependence 
of the dark energy density on the redshift z can be derived from Eq. ([7]) . After a simple algebraic calculation, one 
can obtain 



p = P i 1 ± (1 - Ay + 6 7 ln(l + z) 



1/2 



, A = po/pt, 7 = a 2 /x 2 . 



(54) 



The sign "+" corresponds to the case of quintessence (A > 1) while sign "— " to that of phantom energy (A < 1) 
The current EoS parameter wq is 



w = -1 - 



1 



A(l-A)" 
Therefore, the dependence of the luminosity distance Dl on the redshift z is 



(55) 



d l = -£-(i + z) / (n m (i + zf + n B h{z)) 1/2 dz, 

H o Jo 

h(z) = A" 1 ^1± ((l-A) 2 + 6 7 ln(l + z)) 1/2 ") . 



(56) 
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Table II: Numerical estimation of the difference tt — to in Gyr for various values of A and 7 in the case of a phantom model. 



A 


TOO = -1.10 


too = -1.08 


too = -1.06 


too = -1.04 


too = -1.02 


0.5 


22.6 


28.4 


38.2 


57.6 


116.1 


0.75 


7.3 


9.9 


12.5 


20.3 


39.8 


0.95 


1.1 


1.5 


2.0 


3.1 


6.3 



Table III: Numerical estimation of the difference ti — to in Gyr for various values of A and 7 in the case of quintessence. 





too = -0.98 


too = -0.96 


too = -0.94 


too = -0.92 


too = 


= -0.90 


A 


At d 


Ai t 


At d 


At f 


Af d 


Af f 


At d 


Atf 


At d 


Ai t 


1.05 


5.78 


5.79 


2.8 


2.81 


1.84 


1.86 


1.36 


1.38 


1.08 


1.10 


1.25 


26.70 


26.74 


12.88 


12.96 


8.34 


8.45 


6.10 


6.24 


4.76 


4.94 


1.5 


46.32 


46.42 


22.52 


22.72 


14.57 


14.85 


10.60 


10.96 


8.23 


8.65 



Eq. (|56p coincides with (jT4j if 7 = (f(x) = 0). The SNe data are available in the range < z < 1.5. Therefore, if 
the parameters A and 7 are such that (1 — A) 2 <C 67111(1 + z) in the observable range, the model under discussion is 
indistinguishable from ACDM cosmology. 

The time remaining before a future singularity is 



if 



1 f" 
f o = 77" / M~L + u)- 1 (fi m (l + uf + Qvhfu)) 

n Ju 



-1/2 



(57) 



The variable u = ao/a — 1 varies from (present time) to exp(— (1 — A) 2 /67) — 1 (p — > 00). The function h(u) coincides 
with h(z) in Eq. (|56j) (after changing z — > u). Numerical estimation of the difference if — io for various values of A 
and 7 is given in Table [TT] and Table HTfl For quintessence, the difference id — io (id is the moment of time when a = 
and deceleration begins) is also calculated. We use the value of the Hubble parameter Hq 1 = 13.6 Gyr. 

The difference 8^1 = 5\og(D/D sc ) (/i is the distance modulus) for < z < 1.5 lies in the interval (—0.35 to 0.35) 
for these parameter values. Taking into account that errors in the definition of the SNe modulus arc ~ 0.075 — 0.5, 
we conclude that our model fits these data with excellent precision. Note that futher observational support for 
quintessence which leads to a Type II future singularity is given in Ref. [24 1 . 

The deceleration and jerk parameters are given by 



<?o = 
Jo = 



97 3 

2A(1-A)" D + 2 a 



1 



97 



2A(1-A) 



1 



7 



(1-A)* 



n D + i- 



(58) 



It is convenient to present qo and jo through the parameters wq and A: 

Qo 



9 3 



Jo 



2 {W ° 



1)1- 



A 



1 



■(tuo + 1) Od + 1 



(59) 



For too ~ — l, the jerk parameter jo differs significantly from the standard value only for A — > l. The deceleration 
parameter go for —1-05 < w < —0.95 lies in the interval qy C — 0.16 < qo < qo C + 0.16. Taking into account the errors 
in the definition of go and jo one can conclude that dark energy models with a possible big crush or sudden future 
singularity fit well the current observational data. 

The example considered above is a good theoretical illustration of dark energy models mimicking vacuum energy 
but leading to singularities of Type II. Dark energy with such behavior can be realized if the function f(x) has a 
singularity at x = Xf . 

An important remark is in order. For quintessence the disintegration of bound structures before a future singularity 
seems to be impossible. From Eqs. (|5|) and ([9]) it follows that 



-Finer = —ml— \ W -\ 

2 I -r 3 



(60) 
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Therefore, the maximal value of the inertial force for quintessence is 

-max = ™±P 
iner o 

The energy density of quintessence decreases when the universe expands and the inertial force also decreases with 
time. If we consider an EoS for which f(p) changes sign at x — x p h then the energy density increases with time and 
tends to x^ h . Hence, such an expanding universe is a de Sitter one. 

VII. BIG CRUSH DARK ENERGY MODELS 

Another way to construct cosmological models with various types of evolution is to define the Hubble parameter 
as a function of time. Let us consider the Type II singularity, for example. This case occurs when H is finite and H 
diverges but is negative. In this case, even though the universe is expanding, all structures are crushed rather than 
ripped. An example is given by 

H = # (0) + H[ 0} (t c - t) a . (61) 

Here Hq and H{ arc positive constants and a is a constant with < a < 1. If we choose a to be given by the 
inverse of an odd number, a = l/(2?i + 1), with positive integer n, we can extend H beyond t = t c by defining 



H^ + h[ 0) (t c -t) a when t < t c 
H (o) _ H (o) q _ t ^ a when t > u 



H - \ rr(0) rr(0) ,, , ,a , , . , " ( b2 ) 



In a sense, the H obtained here is a smooth function of t although H diverges at t — t c since there is no sharp point, 
that is, a point where the line folds. At the point H diverges but H is the continuous, single valued, and monotonically 
increasing or decreasing function of t. One may consider the following model 

J H^ ] + H[ 0) (tanh ^l) a when t < t c 
H= { i4 0) - < J (tanh ££1) a when t > t c ' ^ 

In the model (f6"5)l . we find H — > const when t — > ±oo, that is, the space-time is asymptotically de Sitter. Note that 
when H is finite, a rip occurs only for the phantom case, since H > before the singularity. However, a crush occurs 
for quintessence with H < 0. In the models (|62[) and (|63[) . a rip occurs when H[ < and a crush occurs when 

^ 0) > 0. 

From the FRW equations Q , the total energy density and pressure of dark energy as functions of time are: 

p(t) + p m (t) = 3H 2 , P (t) = -2.H 2 - 2H . (64) 

Neglecting the matter density, one can easily obtain the EoS of dark energy corresponding to model (|52"j) : 

_ / -p + 2a3 n H[ Q) (p 1 / 2 - p f 1/2 )- 2 " , pi = 3i4 0)2 when t < t c 
P \ -p 2a3" J ff 1 (0) (p 1 /2 _ p J/a ) -2„ t pi = 3^(0)2 when t>tc ■ m 

When H^ < or H{ > 0, p — > — oo (phantom Type II singularity) or p — > +oo (quintessence Type II singularity) 
at t — > t c — 0. The scale factor is finite at t — > t — 0. 

a(t) = ao cxp{^ 0) t + ^ 0) (t c -<) Q+1 /( a + l)-^ 0) C +1 /( a + l)} . (66) 

In the following, we assume a is given by a = l/(2n + 1) with positive integer n. We now investigate if any object 
can be ripped or crushed at the singularity t = t c although the inertial force -Fi ner diverges. For this purpose, we 
consider the work or the shift of the kinetic energy of the particle. For purposes of this estimation, we neglect all 
forces aside from the inertial force, and we neglect the first term in ([5]) since we assume only H diverges. Then by 
solving the equation of motion 

rnx = -Finer ~ rnlH = —mlH\ a (t c — t) a , (67) 
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for the model (|62|) . one finds 

lH (o) 

x = x + v t 1 (t c -t) a+1 . (68) 

a + 1 

Then the shift of the kinetic energy can be estimated to be 

AT = J Fidt ~ -mv H[ 0) a f (t c - tf' 1 dt + ml 2 H[ 0)2 a f (t c - t) 2 ^ 1 dt . (69) 

Since a > 0, the integration and therefore AT is finite. Hence if the magnitude of the binding energy or the energy 
supporting the bound object is larger than the absolute value of AT. the object is not ripped or crushed although 
T; nor is infinite at t = t c . 

In the case of a Type III singularity, H behaves as in (|6"Tj) but a is negative and greater than unity, — 1 < a < 0. 
We also note that H [ ' > since H > when t < t c . Then in the inertial force ©, the first term behaves as 
H ~ (t c — t) a ~ and the second one as H 2 ~ (t c — t) a . Since — 1 < a < 0, the first term dominates. Then in a way 



similar to (|69|) . the shift of the kinetic energy can be estimated as 

AT ~ -mv H[ 0) (t c t) a + m< 7 (t c tf a , (70) 



;2 W (0) 2 



when t < t c and it becomes positive and diverges when t — > t c . Therefore the Rip surely occurs even for a Type III 
singularity, which is different from the Type II singularity in (|69|) . where all the objects are not always crushed or 
ripped. 

By using the formulation in Rcf . [25J , we now consider what kind of scalar tensor model, whose action is given by 

dix ^9 [^R- 1^)8^8^ - V(cf>)\ , (71) 

can realize the evolution of H given by Eq. (fBTj) . Here, uj(4>) and V(<j>) are functions of the scalar field <f>. If we consider 
the model where u>(4>) and V(<fi) are given by a single function /(</>), as follows, 

«(0 = -4/"(0) , V{<t>) = ^ (3/W + f "{<!>)) , (72) 

the exact solution of the FRW equations has the following form: 

4> = t, H=f{t). (73) 

Then for the model fl61]) with a = l/(2n + 1), we find 

W (0 = 2 JL-^ {tc - 4>Y^ , V(cf > ) = ±{(H { a) +H[ Q \t c - ( f > y^y-2H[ 0) a(tc-^-^ . (74) 

If we redefine the scalar field as 



(2n + l)i?j 0) 



V = - V , ^ n (*c ~ &^ , (75) 

«;(n + 1) 



the kinetic term in the action (|71|) becomes canonical 



and the potential is given by 




-iw«>)a„*a»* = -ia„»>c>>>, (76) 



'2if 1 (0) (2n+l) 



H^a 


«(n+l) 


K 2 


t yW^n + l) 


v(4>)~ 


'4 2 ^ ' 
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Note that y> < when <f> = t <t c and <p — > when = i — >• t c . Near the singularity <f> = t —* t c (tp — > 0), only the last 
term in the potential (|77[) dominates: 

F(0) ~ -qJi | - ; ^^ y I . (78) 

In particular, when n = 1, we find 

~ /qet(0) 3 

(79) 

Thus, the big crush occurs when the scalar field drops into the infinitely deep potential proportional to the inverse 
power of the scalar potential. 

We have constructed models which generate the big crush and have given the explicit action in terms of the scalar 
field. After the big crush, the universe may evolve to asymptotic de Sitter space-time. Hence, big crush phenomenon 
looks much less dangerous than disintegration of bound structures. 

VIII. PHANTOM MODELS AND SINGULARITIES 

One way to realize many of the models proposed here is through a scalar field with a negative kinetic term (phantom 
models). The asymptotic future evolution of such models was examined systematically in Ref. [26j, and we restate a 
number of those results here in order to show explicitly the relation between various types of future singularity. 

The simplest phantom models are characterized by a field with a negative kinetic term. Such models evolve 
according to the equation 

4> + 3H<j)-V'(^) = 0, (80) 

where the prime denotes the derivative with respect to 0. A field evolving according to this equation rolls uphill in 
the potential. The density and pressure for the phantom field are given by 



¥ 



2 



V{4>), (81) 



and 



respectively, so the equation of state parameter is 



!l2 



-^-V{4>), (82) 



(i/2)? + y(fl, 

(1/2)02 _ v{(j)) 

As noted in Ref. |26| , the asymptotic behavior of the equation of state parameter depends on the corresponding 
asymptotic behavior of V'/V. If V'/V — > 0, then w — > — 1. This set of models displays the most diverse behavior, 



since it can correspond to either a Big Rip, a Little Ri 
V{4>). A Big Rip (Type I singularity) occurs when [26 1 



d, or a Pseudo-Rip, depending on the exact functional form for 



^gW-> finite. (84) 

If, instead, the integral in equation (|8"4"f diverges, we have either a Little Rip or a Pseudo-Rip. A Pseudo-Rip occurs 
if V(4>) — > const, while a Little Rip occurs if V{4>) — > oo (see also Section ITT|) . 

The second set of models examined in Ref. [26j corresponds to V'/V — > constant. This gives a constant value for 
w with w < — 1, and produces a Big Rip (Type I) singularity. 

Finally, if V'/V — > ±oo, we have w — > — oo, which can result in a Type III singularity (see also Ref. [271]). 
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IX. CONCLUSION 

In summary, dark energy models with various scenarios of evolution have been presented. Specifically, we con- 
structed scalar dark energy models with Type II and Type III finite-time future singularities, Little Rip and Pseudo- 
Rip cosmologies with finite-time disintegration of bound structures and Big Crush cosmologies. It was shown that 
such models are consistent with observational data from the Supernova Cosmology Project and therefore may be 
viable alternatives to the ACDM cosmology. Moreover, they may be stable for billions of years before entering a soft 
future singularity (with a finite scale factor at the Rip) or before entering a finite-time dissolution of bound structures. 

We have shown that the future evolution of the universe is determined by the selected EoS of dark energy. Unfortu- 
nately, current data for such important parameters as qo and jo are not very reliable, so the nature of the dark energy 
cannot yet be determined, and one can therefore only consider some typical models. In the future, more accurate 
measurements of the deceleration and jerk parameters as well as other cosmological parameters will help to define the 
exact nature of dark energy. Then we will acquire the information on the parameters of the fluid description for the 
EoS of dark energy in this paper and therefore also the information of the parameters in a reconstructed scalar field 
theory. In other words, more precise values of cosmological parameters may significantly constrain the dark energy 
models under discussion. 

The key point is that the current observational data do not answer, even in principle, the question of whether or 
not the universe will end in a future singularity or Rip cosmology. One can construct (as we have here) models that 
mimic standard ACDM up the present, but evolve in the future into any number of possible future states, including 
Pseudo-Rip models, Little Rip models, and a variety of different future singularities. With a variety of ACDM-like 
cosmological models in hand, one can already start to think about future cosmological experiments to define the 
future of the universe more precisely. 
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